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Two-body decay widths of lowest lying and next-to-lowest lying scalar and pseudoscalar mesons are studied in 
Generalized Linear Sigma Model (GLSM) of low-energy QGD. This model which considers mixing between “two 
quark” and “four quark” chiral nonets has been employed to investigate various decays and scatterings in low 
energy region of QGD. In this paper, r[/o(980) KK] and r[ao(980) —)• KK] are obtained and it is shown that 
two-body decay widths of lowest lying mesons are well predicted by this model while for the next-to-lowest lying 
mesons, only some of the decay widths agree with the experimental results. We have compared the predicted 
decay widths in GLSM with the results obtained in single nonet linear sigma model (SNLSM) to indicate that 
chiral nonet mixing greatly improves the predictions of SNLSM for decay widths. 
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1. Introduction 

Since low energy QCD is nonperturbative, the usual field theory methods are impractical and researchers 
have been examining effective field theory approaches. In theses approaches like chiral perturbation 
theoryPt^ linear and non-linear sigma models^^i^ extended linear sigma modeP^®^ and generalized 
linear sigma model (GLSM) the underlying chiral symmetry of QCD is used in order to construct 
hadronic Lagrangians. 

Generalized linear sigma model is successful in describing some unusual features of scalar and pseu¬ 
doscalar mesons below 2 GeV by allowing mixing between a quark-antiquark nonet and a four-quark 
nonet. As we know, the light scalar nonet has unusual spectroscopy compared to the vector one. The 
vector mass ordering corresponds to the number of s-type quarks in each state, i.e., m[p(776) ^ ud] < 
?7i[Ar*(892) ^ us] < m[(/)(1020) ^ ss], while this is not the case for the scalars. In fact, the conventional qq, 
SU (3) picture doesn’t match with the mass of the scalar mesons below I GeV. For example, the isodou¬ 
blet ArQ(800) should be heavier than the isovector ao(980) but: to[A'o(800) ^ us] < m[ao(980) ^ ud]. 
In 1977, Jaffe examined the possibility that the light scalar mesons are made of two quarks and two 
anti-quarks (MIT bag model) and showed that this proposal can explain the unusual mass spectrum of 
the scalar mesons.l^iH^ In this picture the mass ordering of the scalar mesons becomes proportional to 
the number of s-type quarks in each state. As an example, Kq with one strange quark becomes lighter 
than ao(980) with two strange quarks, m[A'Q(800) ^ dsud\ < m[ao(980) ^ dssu], which is consistent with 
the experimental data. Others also studied properties of light scalar mesons in tetraquark picture (see 

Refs.lMMI- 

Since this pure four-quark picture is deficient in describing some decay properties of lowest lying 
scalars and besides the fact that the pure quark-antiquark picture for the next-to-lowest lying scalars can 
not completely describe their properties, mixing among these states was considered as a solution to this 
puzzling feature.l^^lt^ Mechanism for such mixing in a nonlinear chiral Lagrangian was studied in Ref. 
[33 and it was shown that due to level repulsion, the lowest lying scalar mesons become lighter and this 
gives an answer to the question of why scalar mesons are lighter than the axial vector mesons. 

The effect of possible mixing among two quark and four quark chiral nonets on the masses and some 
decay properties of scalar and pseudoscalar mesons, was studied in generalized linear sigma model and it 
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was found that the scalar mesons below 1 GeV are predominantly four-quark states while the lowest lying 
pseudoscalar mesons are closer to qq typeP^ The role of these mixing patterns were also explored in tttt 

and scatterings and it was shown that the model prediction for the scattering amplitudes 
are reasonable and in good agreement with the experiment for the energy region below 1 GeV. Also the 
poles in the scattering unitarized amplitudes, which correspond to the physical resonances, were found 
to represent the scalars with masses and decay widths close to experimental data. Furthermore, the 
r]' —^ rj'K'K decay was probed within this model and it was shown that the predicted decay width agrees 
with the experiment up to about 1 %.!^ 

In this article we will investigate two-body decays of lowest lying and next-to-lowest lying scalar and 
pseudoscalar mesons. For more convenience the experimental data for masses and decay widths are given 
in Table [Tpl A brief review of the model is given in Sec. H] and in Sec. 12] the predictions of the model 
for all hadronic two-body decay widths of scalars and pseudoscalars below 2 GeV are presented. Sec. |4| 
is devoted to the unitarity corrections and their effects on decay widths. In Sec. |5l we summarize and 
discuss the results. 


2. Brief review of the generalized linear sigma model 

The effective Lagrangian of GLSM is constructed out of 3x3 matrices M and M'PSl 


S'! + i(j)\ 

si + i(f>l 

si + i(t>\ 


-k icji'i 

S'l + 


S 2 +m 

Si + i4>l 

S 2 + 

, M' = s'+ icj)' = 

S '2 + 1^2 

s'i + 

S '2 + 


Sl + i4>l 

S'! -k 


S'^ + 

s'i + i4i 

s'i + 


( 1 ) 

where M contains 18 “bare” quark-antiquark scalar and pseudoscalar fields and M' includes 18 “bare” 
scalar and pseudoscalar fields containing two quarks and two antiquarks. Although under chiral transfor¬ 
mations SU(3 )lX SU(3)r, M and M' transform in the same way 

M —^ UlMU^, M' —^ UlM'U^, (2) 

they transform differently under U(l)^ 

M —^ M' —^ (3) 


The Lagrangian density of the model which respects chiral symmetry SU(3 )l x SU(3)R(but not necessarily 
U(1 )a symmetry), is 

£ = -^Tr - iTr - Vo (M, M '), (4) 

where Vo{M,M') at the leading order N < 8 (eight or fewer underlying quark plus antiquark lines at 
each effective vertex) reads 


Vo = - C 2 Tr(MM+) -k 4 Tr(MM^MM+) 

+ d 2 Tr(M'M'+) + + H.c.) 


C3 


7 iln( 


detM 

detAff 


) + (1 - 7i)ln 


Tr(MAL't) 


Tr(Af'Aft) 


(5) 


All the terms except the last two are invariant under U(1)a. We have omitted a possible term [Tr(AfAfl)] ^ 
because it violates the Okubo-Zweig-Iizuka rule. Also we should add a simple chiral symmetry breaking 
term due to the small light quark masses to the potential Vo 


VsB = -Tt[A{M + Aft)] = -2Tr[AS], A = diag(Ai, A 2 , A 3 ), ( 6 ) 


where Ai, A 2 and A 3 are proportional to the three light quark masses. The equilibrium point (or ground 
state) of the system can be found by imposing extremum conditions 
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Table 1. 


Brief review of lowest lying and next-to-lowest lying scalar and pseudoscalar mesons properties 


nn 


particle mass (MeV) full width (MeV) hadronic two-body branching ratios or fractions Ti/Tj 

decay modes 


below 1 GeV 


ttO 

i-(o-+) 

134.9766 ± 0.0006 

~ 10-® 


i-{o-) 

139.57018 ± 0.00035 

~ 10-14 

K± 

1 

o 

493.677 ± 0.016 

~ 10-14 

RO 

^(0-) 

497.614 ± 0.024 

— 

/o(500) or a 

0+(0++) 

400 to 550 

400 to 700 

V 

0+(0-+) 

547.862 ± 0.018 

(1.31 ±0.05) X 10 

/(■^(SOO) or K 

1(0+) 

682 ± 29 

547 ± 24 

rj' 

+ 

1 

O 

+ 

O 

957.78 ±0.06 

0.198 ±0.009 

ao(980) 

l-{0++) 

980 ± 20 

50 to 100 

/o(980) 

0+(0++) 

990 ± 20 

40 to 100 

above 1 GeV 

r;(1295) 

0+(0-+) 

1294 ± 4 

55 ±5 

7r(1300) 

l-(0-+) 

1300 ± 100 

200 to 600 

/o(1370) 

0+(0++) 

1200 to 1500 

200 to 500 


t;(1405) 

o+(o-+) 

1408.9 ±2.4 

51.1 ± 3.2 

if^(1430) 

1(0+) 

1425 ± 50 

270 ± 80 

ao(1450) 

l-(0++) 

1474 ± 19 

265 ± 13 

r;(1475) 

+ 

1 

O 

-t 

o 

1476 ± 4 

85 ±9 

/o(1500) 

0+(0++) 

1505 ± 6 

109 ±7 


/o(1710) 

0+(0++) 1720 ±6 

135 ± 8 

)?(1760) 

0+(0-+) 1756 ±9 

96 ±70 


(20.66 ±0.08)% 


TTTT (dominant) 


ttK 


f TTT] (dominant) 
KK 


r{KK)/r{7Tr]) = 0.183 ± 0.024 


TTTT (dominant) 
KK 


I ao(980)7r 

r(ao(980)7r)/r((jr?) = 0.48 ± 0.22 
ar) 

pir - 


TTTT 


7r(1300)7r 

' 

— 



KK 


j ao(980)7r 


[ /o(980)r, 


ttK 

(93 ± 10)% 

•nri 

V{nri')/V(nri) = 0.35 ± 0.16 

‘ nr]' 

r(i<';^)/r(7r77) = 0.88 ± 0.23 

KK 


ao(980)7r 


TTTT 

(34.9 ± 2.3)% 

7r(1300)7r 

— 

' VV 

(5.1 ±0.9)% 

W 

(1.9 ±0.8)% 

KK 

(8.6 ± 1.0)% 

TTTT 

r(7r7r)/r(i^l?) = 0.4lt°+l 


r{riri)/V(K'K) = 0.48 ± 0.15 
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m\ /dVsB\ 

ds/, \ dS /, 


= 0 , 


dS'/, 


= 0 , 


where the equilibrium values of fields S", (j)', S and (j) are respectively 


{S'th = SlPa, 


irjo = 0 , 


(^^0 = Slaa 


= 0 . 


( 7 ) 

( 8 ) 


In this paper, we assume isotopic spin symmetry so that Ai =A 2 ^ A 3 , ai = a 2 ^ <33 and Pi = P 2 ^ 
p 3 . The four unknown coefficients in the Lagrangian (c 2 , C 4 , ^ 2 , 63 ) and the six unknown parameters 
(ofi, as. Pi, Ps, Ai, A 3 ) can be determinecP^ by using the four minimum potential conditions (Eq. ([7])) 
together with the following experimental input^^ 

m[ao(980)] = 980 ± 20 MeV, 

TO[ao(1450)] = 1474 ± 19 MeV, 
m[7r(1300)] = 1300 ± 100 MeV, 

= 137 MeV, 

JV = 131 MeV, 

As 

—i =20^30, (9) 

Ai 

where A 3 /A 1 is the ratio of strange to nonstrange quark masses!^ Obviously, there are large uncertainties 
in the values of TO[7r(1300)] and A 3 /A 1 which dominate the uncertainty of predictions. 

Since the remaining parameters C 3 and 71 only affect the isosinglet pseudoscalars, one needs 7 masses 
as inputs. From Table [2 it is clear that there are two rj^s below 1 GeV which are good candidates for rji 
and ri 2 predicted by our model 


[77(547)] = 547.862 ±0.018 MeV, 

m®’'P-[7y'(958)] = 957.78 ± 0.06 MeV, (10) 

and four experimental candidates for the two heavier 77 ’s in our model 

m'^’^P-[77(1295)] = 1294 ±4 MeV, 

^exp. [77(1405)] ^ 1408.9 ±2.4 MeV, 

^exp. [77(1475)] ^ 1476 ±4 MeV, 

TO^’^P- [77(1760)] = 1756 ± 9 MeV. (11) 

This leads to six scenarios to identify the two heavier 77 ’s (773 and 774 ) with the four experimental candidates 
above 1 GeV. Tr (M^) and det are the experimental inputs for the determination of C 3 and 71 

Tr{Mj)=Tr{M,7^^. 

det (M“) = det - (12) 

For each of these six scenarios, two sets of 71 and C 3 are found as a result of the quadratic form of 71 
in Eq. (HID. These twelve possibilities are studied in detail in Refs. [20l and l4T] and it was shown that the 
third scenario with the experimental candidates 77(1295) and 77(1760) and solution I (i.e., scenario 31), 
has the best agreement with experimental mass spectrum of the rj system. The variations of these twelve 
parameters are plotted for different values of ?7r[7r(1300)] and A 3 /A 1 in Ref. and it is shown that the 
variations of parameters are most affected by the uncertainty in ?7r[7r(1300)]. 

After fixing the twelve parameters of the model, we will have the rotation matrices describing the 
underlying mixing among two- and four-quark components for each spin and isospin state 


7r+(137) 


Pi 


A:+(496) 


Pi 




5 


1 ^ 

II 


7r+(1300) 


P'l 


A:'+(1460) 


p'l 
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a(J-(980) 


Si 


7fo(800) 


Si 

( 13 ) 

= L-" 


5 

= L-" 


aj(1450) 


S'l 


77^(1430) 


S'l 



where ^ and i?^^are the rotation matrices for 7 = 1 and 7=1/2 pseudoscalars and ^ and ^ are 
the rotation matrices for 7 = 1 and 7=1/2 scalars. For isosinglet scalars and pseduscalars 


h 


fa 

f2 

= Lo^ 

ft 

h 


fc 

_h_ 


fd 


~ 


■ 

m 


Va 

m 


Vb 

V3 


Vc 

V4 


m 


(14) 


where fi, i 
GeV and 


I ..4 and r]i,i = 

...4 are four of the physical 

isosinglet scalars and pseudoscalars below 2 





V2 

cx nn, 

p 

II 

cx nn, 

ft = si 
\ 0^1 1 0^2 

“ V2 

oc ss, 

Vb = (t>l 

cx S5, , , 

(15) 

oc nsnSj 


CX nsns^ 

[fd = S'i 

cx nnnn^ 

Jld = (t>'l 

cx nnnn. 


In this paper, we compute two body decay widths of the scalar and psuduescalar mesons using the 
same order of potential in Ref. [20] with fixed parameters. This provides further test of the underlying 
two and four-quark mixing among the scalar and pseudoscalar mesons below and above 1 GeV and the 
appropriateness of the generalized linear sigma model developed in Ref. |20| and references therein. 


3. Two-body decays 

In this section, we present the prediction of GLSM for hadronic two-body decay widths of lowest lying and 
next-to-lowest lying scalar and pseudoscalar mesons to show whether the mixing of scalar and pseudoscalar 
mesons can improve the results of SNLSIVP^ for the lowest lying decays. We also report the results for 
the decays above I GeV which are beyond the prediction of single nonet model. The decay widths are 
obtained through the following formulas 



where 9 , the center of mass momentum, is given as 


q = - (m 2 -I- m^YWml - {m2 - m.3)^], 


2mi 

for a general two-body decay 1 —>■ 2 3 and the coupling constants are defined as 

n ^ J* I ^fiKK p I I , TJ \ 

- £ = ^ ^ 

+ lKiK'n [kiKt] -I- H.c.) -I- jKiKr,' {kiKt]' -|- H.c.) -f jaj-KT,aj ■ Trry -I- ' W 

+ ihmh'n'n + 7/iw'/iW + + ■ ■ ■ , 


(17) 


(18) 
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where the subscripts i{= 1, 2, 3 and 4), j{= 1, 2) and l{= 1, 2) show the different isosingle, isovector and 
isodoublet meson states, respectively. These isomultiplets contain the physical fields 



K+ 


K-lf 


K,+ 

K = 


K = 

, K = 











1 z 

TTi = — (7r++7r"), 7r2 =-^(7r+- TT”), 7r3=7r°, 

1 z 

Ooi = -^(0+ +Oo), 002 =- 0(7), 003 = o[;. 

The coupling constants are related to the bare couplings through the following relations 


T/i'TTTT — 


1 / d^V 


1 


\ dfi dn+ dn- / ^2 \ dfi 5(0^)^ 


E 




1 / 93 ^ \ 1 

tt' - ^ \ TT !! ^ : 


E 


d^V 


d^V 


IfiKK 


\ dfi dTT+ dn- j V2 \ dfi d{(l)l)A d{(l)l)B 

dfid{(t>l)Ad{(l>\)B 


= V2 


d^V 


dfi dK+ dK- 
d^V 


3 =V^E 


I.A.B 


d^V 


(Eo)/z (-^7r)Al (-^7r)si; 


(Eo )/2 (^7r)Al (-^7r)s2; 


(Eo)/z (-Rk)bi, 


^ajTZT} — 


Jaj nr}' 


^KiKn — 


laiKK = 


d^V 


^ da^ dTT+ drj j \ d{Sf)A d{4>l)B drji 

^ daj dTT+ drj j \ d{Sf)A dr]i 

dnf dK- d7T+ / “ \ d{(t>\)B d{4>l)c 


d^V 


da+ dK° dK- 


= E 


d^V 


AM,C A{Sl)Ad{4>l)Bdmc 


{La)Aj {Rif) Bl (-Ro)n, 

{La)Aj {R-if) Bl {Ro)i2, 

{Lh)aI {Rk)b1 {Rtv) Cl, 
{La)Aj {Rk)bi {Rk)ci, 


^ 1 / d^V \ ly^/_ ^ _ 

hrir, 2 \ Of, dr]dr]/ 2 \ dfidrjjdijK 


{Lo)ii (i?o) Ji iRo)Ki, 


d^V \ ^ / d^V 

IfiW — \ A A 


dfi dr] dr]' / \ dfidr]jdr]K 


{Lo)ii{Ro)ji{Ro)K2- 


(19) 


( 20 ) 


where A, i?, C = 1, 2 with 1 denoting nonet M and 2 denoting nonet M' and I, J, K = o, 6, c and d 
represent the four bases in Eq. (Ha. To, i?T, La, .Ro, -bre, and Rk are the rotation matrices defined in 
Eqs. (11311 and (|14ll . The nonvanishing bare three-point coupling constants are given in Appendix A. The 
contour plots of all the hadronic two-body decay widths below 2 GeV are presented in Figs. [Hill In the 
following we will explain each figure separately: 


Figure [T] 

In this figure we review the results for a -A tttt, /o(980) — >■ tttt, ao(980) —>■ nt] and k -a ttK decay 
widths which are given in Ref. 1411 Contour plot for a -A tttt decay width shows that for a large part of 
the parameter space, the prediction of the model for this decay overlaps with the experimental data (The 
region inside the dashed curves show the experimental range for each decay). Also for fg (980) —?> tttt and 
ao(980) -A nr] decay widths, there exist regions of parameter space which agree with the experimental 
ranges. No region coincides with the experimental range for k -a nK decay width, but for TO[7r(1300)] 
near 1.4 GeV, decay width reaches 400 MeV which at least has a right order of magnitude compared 
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r(a —>7T7T) 



0.693 

0.630 

0.567 

0.504 

0.441 

0.378 


r(^,(980)^7T7T) 
7. 



1.20 1.25 1.30 1.35 1.40 

m[7T(1300)] (GeV) 


0.078 

0.065 

0.052 

0.039 

0.026 

0.013 


r(/f^7TK) 



r(ao(980)^7T77) 



Fig. 1. Contour plots of the prediction of the model for r[cr — > tttt] (top left), r[/o(980) —>■ Trvr] (top right), r[K —>■ ttK] 
(bottom left) and r[ao(980) —^ 7777 ] (bottom right) over the m[7r(1300)]-j43/j4i plane. The parameter spaces inside the 
dashed curves indicate regions for which decay widths are in the experimental range. 


to PDG data. The interesting point is that due to final state interactions in -kK scattering, total decay 
width for k, gets larger (It is shown in Ref. 131] and will be discussed in the next section). 


Figure [2] 


Since tttt and irr] decay modes are dominant channels for decays of /o(980) and ao(980), respectively, 
KK decay channel for these two particles should have narrow decay width. Fig. j^jshows that the model 
prediction for these channels fulfill our expectation: decay widths are small for some regions of the 
parameter space. Averaging over the entire parameter space, we have 

r[ao(980) ^ KK] = 21 ± 27, 

r[/o(980) ^ ATAT] = 21±31, (21) 


where the uncertainty is the standard deviation of the predicted data. These values are in agreement 
with the results of Refs. HSlandHH The prediction of the model over the TO[7r(1300)] — A 3 /A 1 plane for 
r[ao(980) ^ A:i^]/r[oo(980) ^ nri] is 


r[ao(980) ^ KK] 
r[ao(980) —>■ 7r?7] 


0.44 ±0.41. 


( 22 ) 


The mean value is about 2.5 times larger than the experimental value (Table |T]) but the uncertainty 
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r(ao(980)^KA') r(/&(980)^KA') 

30 

28 

26 


22 

20 

1.20 1.25 1.30 1.35 1.40 1.20 1.25 1.30 1.35 1.40 

m[TT(1300)] (GeV) m[7T{1300)] (GeV) 

Fig. 2. Contour plots of the prediction of the model for r[ao(980) —>■ KK] (left) and r[/o(980) —> KK] (right) over the 
m[7r(1300)]-j43/Ai plane. The white region belongs to zero value of decay width. As expected, decay widths are small for 
the main part of the parameter space. 



covers the observed range. 


Figure [3] 

Now let us study hadronic decay widths above 1 GeV. It is shown in Ref. [^Hlthat the GLSM prediction 
for the third f-meson mass is 1504±6, but due to unitarity corrections in tttt scattering,!^ its mass reduces 
to 1149 ± 43 which is close to /o(1370) mass. Furthermore, we know that broad states receive more con¬ 
tributions (in mass and width) from the unitarity corrections compared to the narrow ones. Therefore we 
identify the third f-meson predicted by our model with the broad state /o(1370) and not the narrow state 
/o(1500). This is consistent with the result of Ref. 071 in which it was shown that /o(1500) is the scalar 
glueball state. Moreover, /o(1370) was shown to be predominantly a quark-antiquark state (particularly 
ss) with small remnant of four quark and glue components.!^ Although there is no experimental data 
for decay widths of different modes of this broad resonance. Fig. [3] shows that for all the decay modes 
except TTTT, the predicted decay widths are reasonable compared to the full width but r[/o(1370) — tttt] is 
extremely large and far from experiment. The predicted averaged value for r[/o(1370) — KK] is 72 ± 76 
MeV in agreement with Ref. l48l 


Figure [4] 


For /o(1710) which decays predominantly into kaons, KK and 7777 channels have large and unreasonable 
widths but their ratio is 


r[/o(i7io) ^ 7777 ] 
r[/o(1710) ^ KK] 


0.42 ±0.04, 


(23) 


which overlaps with the experimental range (Table [T]). For tttt channel, decay width has the right order 
of magnitude compared to the full width for large values of 777[7r(1300)] and 


r[/o(i7io) ^ H 
r[/o(1710) ^ KK] 


0.10 ±0.05, 


which is not in the range of PDG data but consistent with the results of Ref. 051 


(24) 


r[/o(1710) ^ 
r[fo{1710)^ KK] 

Moreover the recent results of Refs. Eni and m have suggested that the /o(1710) is predominantly the 
glueball state and due to the fact that the glueball effects are not included in GLSM, we do not expect 
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to get good results for decay widths of this state. 


CO 




r(/^0(1370)^7T7T) 



U hi . . . , I . , . 1 

1.20 1.25 1.30 1.35 1.40 


m[ 7 T( 1300 )] (GeV) 


r(/&(1370)^KA') 



1.20 1.25 1.30 1.35 1.40 

m[TT(1300)] (GeV) 


r(«(1370)^A7A7) 

0.312 
0.260 
0.208 
0.156 
0.104 
0.052 


1.20 1.25 1.30 1.35 1.40 

m[7T(1300)] (GeV) 



CO 




r(^,(1370)^7T(1300)7T) 



0.132 

0.110 

0.088 

0.066 

0.044 

0.022 


Fig. 3. Contour plots of the prediction of the model for r[/o(1370) tttt, KK^ 'qr}, 7r(1300)7r] over the m[7r(1300)]-A3/Ai 
plane. For all decay modes except tttt, predicted decay widths are reasonable compared to the full width. 


Figure [5] 

r[irQ(1430) —>■ ttK] is too large, but according to final state interactions of ttK which are estimated by 
the K-matrix unitarization method, decay width decreases,^ but this time becomes too small compared 
to the experiment (Table 0]) and this convinces us that the current Lagrangian of the model, cannot 
describe energy region above 1 GeV for some processes and therefore is not complete. 


Figure [6] 


Among different hadronic channels for ao(1450) decay, only r[ao(1450) —>■ ttt]'] is comparable with the 
experimental full width and the other two, i.e., irr] and KK widths, are too large but their ratio is in the 
experimental range 


r[ao(1450) ^ KK] 
r[ao(1450) —!> nr]] 


0.96 ±0.10. 


(25) 




















10 S. M. Zebarjad and S. Zarepour 


5 

% 


30 

28 

26 

24 

22 

20 

1.20 1.25 1.30 1.35 1.40 

m[7T(1300)] (GeV) 





30 

0.486 


28 

0.405 



0.324 


26 

0.243 


24 

0.162 



0.081 


22 



20 


U b^ 

1.20 


r(fo(1710)^KA') 



1.25 1.30 1.35 1.40 

m[7T{1300)] (GeV) 


4 

3.5 

3 

2 5 


5 

ro 



1.76 

1.60 

1.44 

1.28 

1.12 

0.96 


Fig. 4. Contour plots of the prediction of the model for r[/o(1710) —>• tttt, KK, r/r)] decay widths over the m[7r(1300)]- 
A 3 /A 1 plane. KK and rjrj channels have large and unreasonable widths, but for tttt channel, decay width has the right order 
of magnitude for large values of m[7r(1300)] compared to the full width. 


r(/<'o‘(1430)^7TK) 



ksi 


5.25 

4.75 

4.25 

3.75 


Fig. 5. Contour plot of the prediction of the model for Kq ttK decay width over the m[7r(1300)]-A3/Ai plane. Decay 
width is too large but receives unitarity corrections due to the ttK final state interaction. 
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Fig. 6. Contour plot of the prediction of the model for r[ao(1450) —y 7 rr}, KK, Trr]'] decay widths over the m[7r(1300)]- 
A 3 /A 1 plane. r[ao(1450) —>■ ttt)'] is comparable with the experimental full width and the other two decay widths are too 
large but their ratio is in the experimental range. 


It should be pointed that the final state interaction effects are also computed for ttt] scatterin^^ and 
fortunately the total width for ao(1450) which comes from the imaginary part of the K-matrix unitarized 
amplitude pole, is close to the experimental range (Table 0]). This will be discussed in the next section. 


Figure [3 


Decay widths of different modes of 77(1295) are calculated within scenario 3/ which is the best scenario. 
Decay widths for both channels are consistent with the experimental data for small values of m[7r(1300)] 
and their ratio is 


r[r7(1295) ^ ao(980)7r] 
r[ 77 ( 1295 ) ^ 0-77] 


0.47 ±0.16, 


(26) 


which overlaps with the experimental range. 


Figure [8] 

To evaluate 77(1405) and 77(1475) —>■ ao(980)7r decay widths, we used scenario 51 and 61, respectively. 
Although 77 masses predicted in these scenarios are less than expected values and these scenarios are 
clearly not favoured, but the resulting widths are comparable with the full experimental width (Fig. 
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0 ). Note that in our model, 77(1405) can not decay into fo{980)r] as its predicted mass is less than the 
threshold mass. 


r[/ 7 ( 1295 )^ao( 980 ) 7 T] 



1.20 1.25 1.30 1.35 1.40 

m[7T(1300)] (GeV) 


0.090 

0.075 
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0.015 
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m[7T(1300)] (GeV) 


^[n{^295)^an] 


— 0.204 
0.170 
0.136 
0.102 
0.068 
0.034 


Fig. 7. Contour plot of the prediction of the model for r[r;(1295) —)• ao(980)7r, err)] decay widths over the m[7r(1300)]-A3/44i 
plane. Decay widths are consistent with the experimental data for smaller values of m[7r(1300)]. 


r[/ 7 ( 1405 )^ao( 980 ) 7 T] r[/ 7 ( 1475 )^ao( 980 ) 7 T] 

0.18 
0.15 
0.12 
0.09 
0.06 
0.03 


1.20 1.25 1.30 1.35 1.40 1.20 1.25 1.30 1.35 1.40 
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Fig. 8. Contour plots of the prediction of the model for r[77(1405) —)• ao(980)7r] and r[77(1475) —)• ao(980)7r] decay widths 
over the m[7r(1300)]-A3/Ai plane. Decay widths are comparable with the full experimental widths. 


In order to compare the predictions of the model with experimental data, the averaged values and 
the standard deviations over the 777[7r(1300)] — A 3 /A 1 plane for decay widths of lowest lying scalars and 
pseudoscalars are given in Table [2] Also the mean values for decay widths of next-to-lowest lying mesons 
are summarized in Table [S] Among these decays, nine of them overlap with the experimental range or 
at least have the right order of magnitude compared to the experimental full width. For the remaining 
decays, i.e., /o(1370) —?> tttt, /o(1710) —>■ KK, /o(1710) — ?> 7777 , ArQ(1430) —>■ ttAT, ao(1450) —>■ 7777 and 

ao(1450) —7> KK, decay widths are too large. 

One may argue that if we consider the effect of final state interactions in tttt, ttK and 7777 scatterings 
within the GLSM framework, this may change the bare predictions of decay widths and masses of scalars 
to the acceptable ones. In the next section, we will review this effect. 
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Table 2. Predicted bare masses and decay widths of low¬ 
est lying mesons in the generalized linear sigma model. 



GLSM 

Width (MeV) Mass of decaying 

particle (MeV) 

Cr — >■ TTTT 

531 ± 99 

645 ± 42 

/o(980) —>■ TTTT 

35 ± 27 

1131 ±47 

/o(980) ^ KK 

21 ± 31 

1131 ±47 

K. —>• TtK 

58 ±90 

1105 ±28 

ao(980) —>■ nr] 

57 ±44 

980 ± 20^ 

ao(980) KK 

21 ± 27 

980 ± 20 


^ This shows experimental error bar; Other error bars 
are standard deviations of the mean value averaged 
over the m[7r(1300)] — A 3 /A 1 plane. 


Table 3. Predicted bare masses and decay widths of next-to-lowest lying mesons in GLSM. 


decay modes 

GLSM 

width (MeV) mass of decaying 

particle(MeV) 

decay modes 

GLSM 

width (MeV) mass of decaying 

particle(MeV) 

/0(1370) ^ TTTT 

3216 ±313 

1504 ± 6 

ao(1450) —r nr] 

1729 ± 278 

1474 ± 19" 

fo{l370)^KK 

72 ± 76 

1504 ± 6 

ao(1450) —>• nr]' 

106 ± 57 

1474 ± 19 

/o(1370) ^ r?r, 

107 ± 98 

1504 ± 6 

ao(1450) —>■ KK 

1660 ± 278 

1474 ± 19 

/o(1370) 7r(1300)7r 

59 ±52 

1504 ± 6 

??(1295) -J- no (980)71 

62 ±30 

1271± 61 

/O (1710) ^ TTTT 

339 ± 204 

1684 ± 40 

»?(1295) <7rj'° 

103 ± 72 

1271± 61 

/o(1710) ^ KH 

3306 ± 523 

1684 ± 40 

77(1405) ^ ao (980)77'= 

65 ±35 

1274 ± 59 

fo(inO)^r]V 

/S'*(1430) nK 

1406 ± 337 

4674 ± 562 

1684 ± 40 

1555 ± 36 

77 ( 1475 ) -7 ao (980)77 

69 ±39 

1278 ± 59 


^ This shows experimental error bar; Other error bars are standard deviations of the mean value averaged over the 
m[7r(1300)] — A 3 /A 1 plane. 

^ Scenario 31 
^ Scenario 5/ 

^ Scenario 67 


4. Unitarity Corrections 

In order to consider final state interactions in tttt, ttK and ttt] scatterings, we use K-matrix unitarization 
methocP^l®l through which the partial wave bare amplitude T/ ^ transforms to unitarized amplitude T/ 
using the following equation 


A 


7^7 B 

AJ 


1 - iT{ 


I 


(27) 


where I and I denote to the partial wave isospin and angular momentum. The poles of the K-matrix 
unitarized amplitude are used to calculate the physical masses and full decay widths of the intermediate 
scalar mesons. Here we review the tttt scatterinj^ and the ttK and irr] scatterings are explored in Refs. 

[Miiini 
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The K-matrix unitarized tttt scattering amplitude for the / = J = 0 channel is given by 


rpO _ 

-^0 — 


7^0-^ 

-^0 


l-zT° 


B ’ 


where is the “bare” scattering amplitude calculated from the Lagrangian in Eq. 0 






lUf^-S 


with 


1 / 4m^ 

Ta =- \ I -^ 

647 r V s 


-5 7ff + 4i:7?„ln(l + l4 

pi \ mi. 








(28) 


(29) 


(30) 


where Ptt = y^s — 4m^/2. The scalar-pseudoscalar-pseudoscalar couplings 7 /^ 7771 - are defined in Sec. [31 and 
is the pion four-point coupling constant 

/ d'^V \ / \ 

^ = \dn+dn-dn+dn-/ = ^ 

(31) 

where the sum is over “bare” pions and A, B, - ■ ■ =1,2 with 1 denoting nonet M and 2 denoting nonet 
M' and Rt; is the pion rotation matrix. The real part of the unitarized scattering amplitude, Eq. (EHI), 
is plotted in Fig [31 It is clear that the prediction of the model is in reasonable qualitative agreement 
with the experimental data up to about 1 GeV. The physical pole positions in the unitarized scattering 



Fig. 9. Real part of the unitarized scattering amplitude of tttt scattering for three different choices of m[7r(1300)] and 
A 3 /A 1 = 30. 


amplitude are determined by solving for the complex roots of the denominator of the K-matrix unitarized 
amplitude Eq. (ESI) 

1 - = 0. (32) 

Each complex pole is equivalent to m| —tmTi, where rhi and T^ are the physical mass and width of the i-th 
pole. By solving Eq. dSlD numerically, we will find four solutions for the pole positions, each corresponds 
to the physical mass and total width of fiS. The masses and total widths of the tttt unitarized scattering 
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amplitude are plotted in Figs. Hn] and [Tl] Comparing the bare masses resulted from the Lagrangian 
(Fig. [12) with the physical masses obtained from the unitarized tttt scattering amplitude iFig. (TUI) , shows 
that the unitarization procedure reduces the bare masses of a and /o(980) which are located out of the 
experimental ranges to physical masses within the experimental bounds. It is clear from Figs. |T] and [TT] 
that before and after unitarization there exists regions of parameter space for which decay widths of 
a and /o(980) cover the experimental ranges 0- Furthermore, for /o(1370) and /o(1710), the bare total 
widths (sum of the average values for different decay modes in Table ED, are too large compared to the 
experimental full width which proves the failure of the model for these decay widths. As we expected, 
the unitarization, makes the total widths smaller with respect to the previous values, but unfortunately 
too small compared to the experimental data IFig. Ill|l . The same unitarization procedure can be applied 
for ttK and irr] scatterings which leads to the unitarized masses and full widths for isodoublets (k, Kq) 
and isotriplets (oq (980), ao( 1450)), respectively. 




Fig. 10. Contour plots of the predictions of the GLSM for the physical masses of four f-mesons obtained from the poles 
of the K-matrix unitarized amplitude of tttt scattering over the m[7r(1300)] — A 3 /A 1 plane. The parameter space inside the 
dashed curves indicate regions for which masses are in the experimental range. 


Finally, the mean values for physical masses and decay widths obtained from unitarization procedure 
corresponding to isosinglet, isodoublet and isotriplet scalar states are given in TablejH It is clear that the 
averaged total widths of cr and ao(980) with their uncertainty overlap with the experiment. For /o(980), 

^As TTTT channel is the dominant decay mode for cr and /o(980) decays, total widths resulted from unitarization procedure 
are comparable with r[cr ^ tttt] and r[/o(980) tttt]. 
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Fig. 11. Contour plots of the predictions of the GLSM for the physical widths of four f-mesons obtained from the poles 
of the K-matrix unitarized amplitude of tttt scattering over the m[7r(1300)] — A 3 /A 1 plane. The parameter space inside the 
dashed curves indicate regions for which decay widths are in the experimental range. 


although the mean value of F/^(gso) is not in agreement with the experimental data, it varies from 75 to 
335 MeV which indeed for some regions of parameter space, total width overlaps with PDG data. F^ does 
not agree with the experiment but is close to it. Thus, the poles in the unitarized scattering amplitude 
which represent physical states leads to the acceptable widths and masses for states below 1 GeV, while 
for the physical states above 1 GeV, the obtained total widths are too small compared to the expected 
values except for ao(1450). 

5. Summary and conclusions 

In the present work, we have studied the effect of underlying mixing between qq and qqqq components 
in hadronic two-body decays of scalar and pseudoscalar mesons below 2 GeV within the GLSM. To see 
whether the mixing effect could improve the prediction of SNLSM for lowest lying scalars,!^ a comparison 
is made between predicted decay widths evaluated in these models ( first two columns of Table[S]). At first 
sight, it is clear that the predictions have greatly improved in GLSM. So this model not only determines 
quark structure of mesons but also can improve the predictions for decays and certainly scatterings of 
mesons 

We should note that the masses of a, /o(980) and their mixing angle 0s, the mass of ao(980) in the 
























Two-body decay widths 17 



1.215 

1.188 

1.161 

1.134 

1.107 

1.080 


'^fo{1270) 



1.20 1.25 1.30 1.35 1.40 

m[TT(1300)] (GeV) 





1.20 1.25 1.30 1.35 1.40 

m[7T{1300)] (GeV) 


1.5081 

1.5048 

1.5015 

1.4982 

1.4949 

1.4916 


1.750 

1.725 

1.700 

1.675 

1.650 

1.625 


Fig. 12. Contour plots of the predictions of the GLSM for the bare masses of four f-mesons over the m[7r(1300)] — A^/Ai 
plane. The parameter space inside the dashed curves indicate regions for which decay widths are in the experimental range. 


Table 4. Predicted physical masses and decay widths of lowest lying and nex- 
t-to-lowest lying mesons in GLSM obtained from unitarized amplitudes.^^^Jt^ 



lowest lying 

Width (MeV) Mass(MeV) 


next-to-lowest lying 

Width (MeV) Mass(MeV) 

CT 

385 ± 61 

476 ±4 

/o(1370) 

22 ±21 

1149 ±43 

/o(980) 

207 ± 65 

1053 ± 44 

/o(1710) 

2 ± 1 

1672 ± 37 

K 

689 ± 27 

722 ± 28 

iiTo* (1430) 

5±6 

1114 ± 50 

ao(980) 

60 ±52 

984 ±7 

ao(1450) 

502 ± 88 

1085 ± 43 


isotopic spin invariant limit and therefore the related vertices such as To-tttt and 7 / 27 r 7 r are not predicted 
in SNLSM directly; The isoscalar masses {m[a], m[/o(980)]) and consequently widths in Table [S] comes 
from the best fit for the real part of the I = J = 0, tttt scattering amplitude. Also the isovector mass 
(m[ao(980)]) can not be predicted from the SNLSM and an arbitrary value should be chosen for it; But 
in GLSM, all the masses and widths are direct predictions of the model and this is another advantage of 
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Table 5. Predicted bare masses and decay widths of lowest lying mesons in the non-renormalizable SU(3) single nonet linear 
sigma model (first column), the generalized linear sigma model (second column) and the decoupling limit of GLSM (third column). 



SNLSM 

Width (MeV) Mass of decaying 

particle (MeV) 

GLSM 

Width (MeV) Mass of decaying 

particle (MeV) 

GLSM in decoupling limit 

Width (MeV) Mass of decaying 

particle (MeV) 

Cr —> TTTT 

830 

847 

531 ± 99 

645 ± 42 

3077=- 

960=" 

/o(980) ->• TTTT 

4109 

1300 

35 ±27 

1131 ±47 

b 

1975 ±117 

/o(980) KK 

1525 

1300 

21 ±31 

1131 ±47 

2424 ± 46 

1975 ± 117 

K, —y TtK 

2350 

1300 

58 ±90 

1105 ±28 

1844 ± 30 

1202 ± 22 

ao(980) —>■ 7r»7 

381 

1100 

57 ±44 

980 ± 20'^ 

347 ± 16 

980 ± 20 

ao(980) KK 

221 

1100 

21 ± 27 

980 ± 20 

161 ± 100 

980 ± 20 


^ These values are independent of A%/A\ and therefore have no standard deviations. 

^ In the decoupling limit /o(980) contains strange quarks and thus cannot couple with tt, i.e., 7/27r7r = 0- 

^ This shows experimental error bar; Other error bars are standard deviations of the mean value averaged over the m[7r(1300)] — 
AzjAi plane in case of GLSM or averaged over different values of A 3 /A 1 in case of GLSM in decoupling limit. 


GLSM versus SNLSM. 

We have also evaluated all the lowest lying decay widths in the decoupling limit 0 {d2, eg ^ 0 
and 7 i —> 1) given in the third column of Tabled in order to show that mixing provides much better 
description of decay widths. It is seen that by decoupling the four-quark fields, the results go far from 
the experimental data. 

As a result, it has become clear that the GLSM Lagrangian at the present order of V = 8 with 8 
experimental inputs and without any further tuning is successful in predicting decay widths of lowest lying 
mesons, while the predictions of SNLSM model are far from the expected values. Therefore, inclusion of 
the underlying mixings considerably improves the results for the decay widths. 

We should notice that GLSM can not give acceptable widths for some states above 1 GeV. It would be 
worthwhile to study the effects of adding terms with higher than eight quark and antiquark lines and also 
considering the effects of scalar and pseudoscalar glueballs in potential. Our initial estimate shows that 
adding scalar and pseudoscalar glueballs has more considerable effects than adding terms with higher 
order of N. 
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Further discussion on decoupling limit is given in [Appendix B| 
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Appendix A. Three-point bare couplings 
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Appendix B. Decoupling limit of GLSM 

The four-quark fields are decoupled in the limit (i 2,63 0 and 71 1. To find five unknown parameters 

in this limit ( 02 , 04 ,^ 1 ,^ 3 , 01 , 03 ), we use Wa and A 3 /A 1 as experimental inputs 

ml = -2c2 -I- 404 O 1 , 

TOq = -2c2 -I- I 204 O 1 , 

= 20 ^ 30, 

Ai 

= 2 oi, (B.l) 


besides two minimum equations 


— 2Ai -|- 4c((oi — 2 oi02 — 0, 

—2A3 -|- 4 o“o 3 — 203O2 = 0. 

The experimental input for the determination of the remaining parameter 03 is Tr(M^ 

ml + m ^ = — 4 o 2 -^ -I- 404 O 1 — -I- 404 O 3 . 

' ' Oi O3 

In this limit, the physical vertices are 

Tcttttt — 4c40;i, 


T/otttt — O 5 

IhKK = —4\/2c40i -I- SV^clas, 
IkttK = 4040^, 


lagirri — 


+4^c;cos(9,)», - 
0 ? of 03 


laoKK = 8C4O1 — 4C4O3, 

where 9p is the pseudoscalar mixing angle. 


(B.2) 

(B.3) 


(B.4) 


V 


cos Op — sin Op 


Va 

v' 


sin Op cos Op 


Ob 


(B.5) 
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